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Abstract. 

In models of triggered seismicity and in their inversion with empirical data, the de- 
tection threshold md is commonly equated to the magnitude too of the smallest trigger- 
ing earthquake. This unjustified assumption neglects the possibility of shocks below the 
detection threshold triggering observable events. We introduce a formalism that distin- 
guishes between the detection threshold md and the minimum triggering earthquake toq < 
TOrf. By considering the branching structure of one complete cascade of triggered events, 
we derive the apparent branching ratio n a (which is the apparent fraction of aftershocks 
in a given catalog) and the apparent background source S a that are observed when only 
the structure above the detection threshold to^ is known due to the presence of smaller 
undetected events that are capable of triggering larger events. If earthquake triggering 
is controlled in large part by the smallest magnitudes as several recent analyses have shown, 
this implies that previous estimates of the clustering parameters may significantly un- 
derestimate the true values: for instance, an observed fraction of 55% of aftershocks is 
renormalized into a true value of 75% of triggered events. 



1. Introduction 

There are numerous evidences that a seismic event can 
have a significant effect on the pattern of subsequent seis- 
micity, most obvious in aftershocks of large events. More 
recently has emerged an important extension of the con- 
cept of earthquake interactions in the concept of triggered 
seismicity, in which the usual distinction, that foreshocks 
are precursors of larger mainshocks which in turn trig- 
ger smaller aftershocks, becomes blurred: a parsimonious 
and efficient description of seismicity does not seem to re- 
quire the division between foreshocks, mainshocks and af- 
tershocks, which appear indistinguishable from the point 
of view of many of their physical and statistical proper- 
ties [Helmstetter and Sornette, 2003a] . An important logical 
consequence is that cascades of triggered seismicity ( "after- 
shocks," "aftershocks" of "aftershocks," ...) may play an 
important role in the overall seismicity budget [Helmstetter 
and Sornette, 2003b; Felzer et al, 2002]. 

There is thus a growing interest in phenomenological 
models of triggered seismicity, which use the Omori law 
as the best coarse-grained proxy for modeling the com- 
plex and multi-faceted interactions between earthquakes, to- 
gether with the other most solid stylized facts of seismic- 
ity (clustering in space, the Gutenberg-Richter (GR) earth- 
quake size distribution and an aftershock productivity law). 
This class of ETAS (Epidemic- Type Aftershock Sequences) 
models introduced by Ogata [1988] and Kagan and Knopoff 
[1981] offers a parsimonious approach replacing the classi- 
fication of foreshocks, mainshocks and aftershocks by the 
concept of earthquake triggering: earthquakes may trigger 
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other earthquakes through a variety of physical mechanisms 
but the effective laws do not allow the identification of a 
particular mechanism. 

The questions suggested by this approach include: 1) 
what is the fraction of triggered versus uncorrelated earth- 
quakes (which is linked to the problem of clustering)? How 
can one use this modeling approach to forecast future seis- 
micity? What are the limits of predictability and how are 
they sensitive to catalog completeness and type of tectonic 
deformation? In general, to attack any such question, one 
needs in one way or another to estimate some key parame- 
ters of the models of triggered seismicity. 

The royal path is in principle to use the maximum likeli- 
hood method to estimate the parameters of the considered 
model from a catalog of seismicity (time, location and mag- 
nitude) (see for instance Ogata [1988] and Kagan [1991]). 
The calculation of the likelihood function requires evaluat- 
ing the theoretical rate of seismicity at time t induced by all 
past events at times ti < t. The maximization of the like- 
lihood with respect to the parameters of the model, given 
the data, then provides an estimate of the parameters. All 
previous studies have considered that small earthquakes, be- 
low the detection threshold, are negligible. Thus, the rate of 
seismicity is calculated as if triggered only by earthquakes 
above the detection threshold. However, this method is not 
correct because it does not take into account events below 
the detection threshold, which may have an important role 
in the triggering of seismicity. Indeed, small earthquakes 
have a significant contribution in earthquake triggering be- 
cause they are much more numerous than larger earthquakes 
(Helmstetter [2003]; Felzer et al. [2002]; Helmstetter et al. 
[2004]). This can simply be seen from the competition 
between the productivity law ~ l0 aM giving the number 
of events triggered by a mainshock of magnitude M and 
the relative abundance ~ W~ bM of such mainshocks given 
by the Gutenberg-Richter law: the contribution of earth- 
quakes of magnitude M to the overall seismic rate is thus 
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~ io-( fc - Q ' M , which is dominated by small M's for a < b 
[Helmstetter, 2003] or equally contributed by each magni- 
tude class for a — b (Felzer et al. [2002]; Helmstetter et 
al. [2004]). Therefore, one needs to take into account small 
events that are not observed in order to calibrate correctly 
models of seismicity and obtain reliable answers to our ques- 
tions stated above. This is an essential bottleneck for the 
development of earthquake forecasts based on such models. 

The purpose of this note is to present a general theoreti- 
cal treatment of the impact of unobserved seismicity within 
the framework of models of triggered seismicity. We show 
by analyzing the branching structure of a complete cascade 
(cluster) triggered by an independent background event that 
the unobserved seismicity has the effect of decreasing the 
real branching ratio n and of increasing the number of in- 
dependent background events S into apparent quantities n a 
and S a - This bias may be very significant. We therefore 
claim that previous work should be reanalyzed from the new 
perspective of our approach. This leads also to important 
consequences for the methods presently used to forecast fu- 
ture seismicity based only on incomplete catalogs. 

2. The ETAS model and the smallest 
triggering earthquake 

2.1. Definition of the ETAS model 

To make this discussion precise, let us consider the 
epidemic-type aftershock sequence (ETAS) model, in which 
any earthquake may trigger other earthquakes, which in turn 
may trigger more, and so on. Introduced in slightly different 
forms by Kagan and Knopoff [1981] and Ogata [1988], the 
model describes statistically the spatio-temporal clustering 
of seismicity. 

The triggering process may be caused by various mecha- 
nisms that either compete or combine, such as pore-pressure 
changes due to pore-fluid flows coupled with stress varia- 
tions, slow redistribution of stress by aseismic creep, rate- 
and-state dependent friction within faults, coupling between 
the viscoelastic lower crust and the brittle upper crust, 
stress-assisted micro-crack corrosion, and more. The ETAS 
formulation amounts to a two-scale description: these above 
physical processes controlling earthquake interactions enter 
in the determination of effective triggering laws in a first 
step and the overall seismicity is then seen to result from 
the cascade of triggering of events triggering other events 
triggering other events and so on [Helmstetter and Somette, 
2002]. 

The ETAS model consists of three assumed laws about 
the nature of seismicity viewed as a marked point-process. 
We restrict this study to the temporal domain only, sum- 
ming over the whole spatial domain of interest. First, the 
magnitude of any earthquake, regardless of time, location, 
or magnitude of the mother shock, is drawn randomly from 
the exponential Gutenberg-Richter (GR) law. Its normal- 
ized probability density function (pdf) is expressed as 
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where the exponent 6 is typically close to one, and the cut- 
offs mo and m m a X serve to normalize the pdf. The upper 
cut-off m max is introduced to avoid unphysical, infinitely 
large earthquakes. Its value was estimated to be in the 
range 8 — 9.5 [Kagan, 1999]. As the impact of a finite m max 
is quite weak in the calculations below, replacing the abrupt 
cut-off m max by a smooth taper would introduce negligible 
corrections to our results. 

Second, the model assumes that direct aftershocks are dis- 
tributed in time according to the modified "direct" Omori 
law (see Utsu et al. [1995] and references therein). Assuming 
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Third, the number of direct aftershocks of an event of 
magnitude m is assumed to follow the productivity law: 



p{m) = fcl0 a(m m °\mo <m< m max 
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Note that the productivity law (3) is zero below the cut- 
off mo, i.e. earthquakes smaller than mo do not trigger 
other earthquakes, as is typically assumed in studies using 
the ETAS model. The existence of the small magnitude 
cut-off mo is necessary to ensure the convergence of these 
types of models of triggered seismicity (in statistical physics 
of phase transitions and in particle physics, this is called 
an "ultra-violet" cut-off which is often necessary to make 
the theory convergent). In a closely related paper, Somette 
and Werner [2004] showed that the existence of the cut-off 
mo has observable consequences which constrain its physical 
value. 

The key parameter of the ETAS model is defined as the 
number n of direct aftershocks per earthquake, averaged 
over all magnitudes. Here, we must distinguish between the 
two cases a = b and a / b: 
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for the general case a=fcb. The special case a = b gives 
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Three regimes can be distinguished based on the value of 
n. The case n < 1 corresponds to the subcritical, stationary 
regime, where aftershock sequences die out with probability 
one. The case n > 1 describes unbounded, exponentially 
growing seismicity [Helmstetter and Somette, 2002]. In ad- 
dition, the case b < a leads to explosive seismicity with finite 
time singularities [Somette and Helmstetter, 2002] . The crit- 
ical case n = 1 separates the two regimes n < 1 and n > 1. 
Helmstetter and Somette [2003b] showed that the branching 
ratio n is also equal to the fraction of triggered events in a 
seismic catalog. 

The fact that we use the same value for the productiv- 
ity cut-off and the Gutenberg-Richter (GR) cut-off is not 
a restriction as long as the real cut-off for the Gutenberg- 
Richter law is smaller than or equal to the cut-off for the 
productivity law. In that case, truncating the GR law at 
the productivity cut-off just means that all smaller earth- 
quakes, which do not trigger any events, do not participate 
in the cascade of triggered events. This should not be con- 
fused with the standard incorrect procedure in many pre- 
vious studies of triggered seismicity of simply replacing the 
GR and productivity cut-off mo with the detection thresh- 
old md in equations (1) and (3) (see, for example, Ogata 
[1988]; Kagan [1991]; Ogata [1998]; Console et al. [2003]; 
Zhuang et al. [2004]). The assumption that m<j = mo may 
lead to a bias in the estimated parameters. 

Without loss of generality, we consider one independent 
branch (cluster or cascade of aftershocks set off by a back- 
ground event) of the ETAS model. Let thus an independent 
background event of magnitude Mi occur at some origin 
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of time. The mainshock will trigger direct aftershocks ac- 
cording to the productivity law (3). Each of the direct af- 
tershocks will trigger their own aftershocks, which in turn 
produce their own, and so on. Averaged over all magni- 
tudes, an aftershock produces n direct offsprings according 
to (4). Thus, in infinite time, we can write the average of 
the total number N to tai of direct and indirect aftershocks of 
the initial mainshock as an infinite sum over terms of (3) 
multiplied by n to the power of the generation [Helmstetter 
and Sornette, 2003b], which can be expressed for n < 1 as: 



Ntotai = p(Mi) + p(M 1 )n + p(Mi)n 2 + ... 
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However, since we can only detect events above the detection 
threshold m d , the total number of observed aftershocks N obs 
of the sequence is simply N to tai multiplied by the fraction 
of events above the detection threshold, given by 
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according to the GR distribution. The observed number of 
events in the sequence is therefore 
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Equation (8) predicts the average observed number of di- 
rect and indirect aftershocks of a mainshock of magnitude 
Mi > md- Sornette and Werner [2004] showed that mo 
may be estimated using fits of Nobs given by (8) to observed 
aftershock sequences and Bath's law. The essential param- 
eter needed to constrain mo is the branching ratio n. As we 
demonstrate below, typical estimates of n in the literature 
obtained from a catalog neglect undetected seismicity and 
therefore cannot be used directly to constrain mo. 

Naturally, there is no justification for assuming that rrid 
should equal mo, as is done routinely in inversions of cata- 
logs for the parameters of the ETAS model (see, for exam- 
ple, Ogata [1988]; Kagan [1991]; Ogata [1998]; Console et 
al. [2003]; Zhuang et al. [2004]). First, detection thresh- 
olds change over time as instruments and network coverage 
become better, while the physical mechanisms in the Earth 
presumably remain the same. No significant deviation from 
the Gutenberg-Richter distribution or the productivity law 
has been recorded as the detection threshold m d decreased 
over time. Second, studies of earthquake occurrence at small 
magnitude levels below the regional network cut-offs show 
that earthquakes follow the same Gutenberg-Richter law (for 
a recent study of mining-induced seismicity, see, for exam- 
ple, Sellers et al. [2003]), while acoustic emission experi- 
ments have shown the relevance of the Omori law at small 
scales (see, for instance Nechad et al. [2004] and references 
therein). Within the assumption of self-similarity, i.e. a con- 
tinuation of the GR and productivity laws down to a cut-off, 
evidence thus points towards a magnitude of the smallest 
triggering earthquake and a Gutenberg-Richter cut-off that 
lie below the detection threshold and are thus not directly 
observable. 

2.2. Two interpretations of the ETAS model 

The ETAS model may be viewed in two mathematically 
equivalent but interpretionally different ways. In this sec- 
tion, we develop both views to underline that our results ap- 
ply in both cases and to stress the equivalence of these two 
views. The first describes the model as a simple branching 
model without loops: The independent background events, 



due to tectonic loading, may each independently trigger di- 
rect aftershocks, each of which may in turn trigger secondary 
shocks, which in turn may trigger more. Because every 
triggered event, excluding of course the non-triggered back- 
ground events, has exactly one mainshock (mother), but 
the mother may have many direct aftershocks (children), 
the model can be thought of as a simple branching model 
without loops. The background events are assumed to be a 
stationary Poisson process with a constant rate. The rate 
of the aftershocks of a background event is a non-stationary 
Poisson process that is updated every time another after- 
shock occurs until the cascade dies out. The intensity is 
thus conditioned on the specific history of earthquakes. The 
expectation of the conditional intensity is an average over 
an ensemble of histories. The predicted number of after- 
shocks of an independent background event of magnitude 
Mi as in expression (8) is thus averaged over the ensemble 
of possible realizations of the aftershock sequence, and it is 
also averaged over all possible magnitudes of the aftershocks. 
The branching ratio n is therefore an average not only over 
magnitudes but also over an ensemble of realizations of the 
non- stationary Poisson process. The model thus consists of 
statistically independent Poisson clusters of events, which 
are, however, dependent within one cluster. 

The second view of the ETAS model does not allow a 
unique identification of the mother or trigger of an earth- 
quake. Rather, each aftershock was triggered collectively by 
all previous earthquakes, each of which contributes a weight 
determined by the magnitude-dependent productivity law 
p(m) that decays in time according to the Omori law tp(t) 
and in space according to a spatial function R(r), often cho- 
sen to be an exponential or a power law centered on the 
event. The instantaneous conditional intensity rate at some 
time t at location r is given by 
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where the sum runs over all previous events i with mag- 
nitude rrii at time U at location n. Thus the triggering 
contribution of a previous event to a later event at time t is 
given by its own weight (its specific entry in the sum) divided 
by the total seismicity rate, including the background rate. 
A non-zero background rate then contributes evenly to all 
events and corresponds to an omnipresent loading contribu- 
tion. In this way, earthquakes are seen to be the result of all 
previous activity including the background rate. This cor- 
responds to a branching model in which every earthquake 
links to all subsequent earthquakes weighted according to 
the contribution to triggering. A branching ratio can then be 
interpreted as a contribution of a past earthquake to a future 
earthquake, averaged over an ensemble of realizations and 
all magnitudes. This second view becomes the only possible 
one for nonlinear models whose triggering functions depend 
nonlinearly on previous events (see e.g. the recently intro- 
duced multi-fractal earthquake triggering model by Ouillon 
and Sornette [2004] and references therein). 

These two views are equivalent because the linear formu- 
lation of the seismic rate of the ETAS model together with 
the exponential Poisson process ensures that the statisti- 
cal properties of the resulting earthquake catalogs are the 
same. The linear sum over the individual contributions and 
the Poisson process formulation are the key ingredients that 
allow the model to be viewed as a simple branching model. 

This duality of thinking about the ETAS model is re- 
flected in the existence of two simulation codes in the com- 
munity, each inspired by one of the two views. A program 
written by K. Felzer and Y. Gu (personal communication) 
calculates the background events as a stationary Poisson 
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process and then simulates each cascade independently of 
the other branches as a non-stationary process. The sec- 
ond code by Ogata [1998], on the other hand, calculates the 
overall seismicity at each point in time by summing over all 
previous activity. The latter code is significantly slower be- 
cause the independence between cascades is not used, and 
the entire catalog is modeled as the sum of a stationary and 
a non-stationary process. Despite the different approach, 
both resulting earthquake catalogs share the same statisti- 
cal properties and are thus equally acceptable. In the in- 
terpretation of the model as a simple branching model, the 
parameter n defined in (4) would correspond to a branch- 
ing ratio, while the view that aftershocks are triggered col- 
lectively by all previous earthquakes and the background 
rate would interpret n as an average contribution of a single 
earthquake on future earthquakes. The important point is 
that the statistical properties are the same. 

While the simulation or forward problem is straight- 
forward when adopting the view of the ETAS model as a 
branching model with one assigned trigger for any after- 
shock, the inverse problem of reconstructing the branching 
structure from a given catalog can at best be probabilistic. 
Because aftershocks of one mother cannot be distinguished 
from those of another mother except by spatio-temporal dis- 
tance, we have no way of choosing which previous earth- 
quake triggered a particular event, or whether it is a back- 
ground event. Rather, we must resort to calculating the 
probability of an event at time t to be triggered by any 
previous event according to the contribution that the previ- 
ous event has at time t compared to the overall intensity at 
time t. This probability is of course equal to the weight or 
triggering contribution that a previous event has on a sub- 
sequent event when adopting the collective-triggering view. 
However, the interpretation remains different since the prob- 
ability specifies a unique mother in a fraction of many real- 
izations. 

Having determined from catalogs a branching structure 
weighted according to the probability of triggering, one may 
of course choose to always pick as source of an event the most 
probable contributor, be that a previous event or the back- 
ground rate. Another option is to choose randomly accord- 
ing to the probability distribution and thus reconstruct one 
possible branching structure among the ensemble of many 
other possible ones. The latter approach has been used by 
Zhuang et al. [2004] and labeled stochastic reconstruction. 

The key point is that equating the detection threshold 
with the smallest triggering earthquake will most likely lead 
to a bias in the recovered parameters of a maximum likeli- 
hood analysis as performed by Zhuang et al. [2004] and in 
many other studies. Therefore, the weights or probabilities 
of previous events triggering subsequent events were calcu- 
lated from biased parameters. 

In the following, we show that the branching ratio and 
the background source events are significantly biased when 
they are estimated from the apparent branching structure 
observed above the detection threshold m d instead of the 
complete tree structure down to m . We adopt the view of 
the simple branching model to make the derivations more 
illuminating but all results can be reinterpreted as contribu- 
tions in the collective-triggering view. 

3. The Apparent Branching Structure of 

the ETAS Model 

3.1. The apparent branching ratio n a 

Seismic catalogs are usually considered complete above a 
threshold m d , which varies as a function of technology and 
location. For instance, m d w 2 for modern Southern Califor- 
nia catalogs (and for earthquakes not too close in time to a 
large mainshock [Kagan, 2003]). The analysis of the statis- 
tics of the Omori and inverse Omori laws for earthquakes of 



magnitude down to 3 (Helmstetter [2003]; Helmstetter and 
Sornette [2003a]) suggests that mo is smaller than the com- 
pleteness magnitude m d and is thus not directly observable. 
Thus, mo is the size of the smallest triggering earthquake, 
which most likely differs significantly in size from the cur- 
rent detection threshold mj. By considering the branching 
structure of the model, we derive the apparent branching 
ratio and the apparent background source that are found 
if only the observed (detected) part of the ETAS model is 
analyzed. 

As stated above, the ETAS model consists of statistically 
independent Poisson-distributed clusters. Within each clus- 
ter, each shock has exactly one trigger, apart from the initial 
background event (mainshock) that sets off the cascade. We 
restrict this study to the case n < 1 for mathematical conve- 
nience and because this range of branching ratios gives rise 
to statistically stationary seismic sequences. 

Since aftershock clusters are independent of each other, 
averages of one cluster are equal to ensemble averages, as 
nothing but the inherent stochasticity of the model deter- 
mines the properties of the clusters. One cluster consists of 
one independent background event (source) and its direct 
and indirect aftershocks. However, if not all events of the 
sequence are detected, then there will appear to be less di- 
rect (and indirect) aftershocks, i.e. the branching ratio will 
appear different. Furthermore, some observed events will be 
triggered by mother-earthquakes below the detection thresh- 
old, resulting in apparently independent background events. 

We will refer to the observed branching structure above 
the detection threshold as the apparent structure. As a 
useful visualization of the effect of the detection threshold 
on the branching structure, one can think of the branching 
structure as a mountain range that is submerged to some 
fraction of its height in water (see Figure 1) . The height of 
each mountain corresponds to the magnitude of an earth- 
quake. To the left (i.e. backward in time) of each peak 
is one peak (the trigger) and to the right (i.e. forward in 
time) of each peak may be several (the direct aftershocks), 
connected from left to right via a ridge linking trigger to 
offspring. Since the height of the peaks varies (according to 
the GR law), some of the peaks will be below water (below 
the detection threshold). A whole part of the branch may 
be submerged in water (i.e. unobserved) until eventually a 
peak rises above the surface and appears to have no ridge 
connecting it to a previous peak, because the ridge of the 
last observable mountain simply descends into the water. 

This view leads to the conclusion that the average num- 
ber of direct aftershocks that are observed will be less than 
the real branching ratio, since some of the triggered events 
of an observed shock will fall below m d and hence not be 
included in the count. Only the fraction f a t s from equation 
(7) above mj of the total direct aftershocks p(m) will be ob- 
served. Moreover, the pdf P(m|m > mj) of mother events 
conditioned on being larger than mj is zero for m < mj 
and equal to P(m)/f b s for m max > m > m d . We can thus 
define the apparent branching ratio as 

n a = J P(m\m > m d )p(m)f obs dm 

m 

m malc 

= J P(m)p(m)dm (10) 
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for the case a =fc b. The special case a = b gives 
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Using equation (4) and eliminating k, we have n a in terms 
of n: 
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when a — b. 

According to expression (12), n a < n, where the equality 
holds for nid equal to mo- In principle, equation (12) also 
holds for n > 1, but we restrict this study to the stationary 
regime n < 1. Figure 2 shows function of n for a 

range of values of mo for the case a = b. The values of 
mo are mo = md = 3 (solid), mo = (dashed), mo = —5 
(dotted) and mo = — 10 (dash-dotted). We assumed mj = 3 
and m ma x = 8. Specifying mo then fixes the linear slope of 
the dependence between n and n a . Figure 2 demonstrates 
that the apparent (measurable) fraction of aftershocks may 
significantly underestimate the true fraction of aftershocks 
even for m not very small. For example, m = — 5 roughly 
translates a real branching ratio of n = 0.9 into an appar- 
ent branching ratio n a = 0.3. Decreasing a below & places 
more importance on the triggering from small earthquakes 
and therefore strongly amplifies this effect. 

In Figure 3, we plot the ratio n a /n as a function of the 
unknown m . Again, we assume nid = 3, m max = 8, b = 1, 
but now we let a = 0.5 (dash-dotted), a = 0.8 (dashed) and 
a = b = 1.0 (solid). As expected, when mo = m^, the ra- 
tio is one because there is no unobserved seismicity. As mo 
goes to minus infinity, n a approaches zero since almost all 
seismicity occurs below the threshold. We see clearly that 
unobserved seismicity results in a drastic underestimate of 
the fraction of aftershocks. 

Given an estimate of the magnitude of the smallest trig- 
gering earthquake mo (see Sornette and Werner [2004] and 
references therein), one can calculate the true branching ra- 
tio from the apparent branching ratio. In fact, Sornette and 
Werner [2004] obtained four estimates of mo as a function of 
n by comparing the ETAS model prediction of the number 
of observed aftershocks (8) from fits to observed aftershock 
sequences and from the empirical Bath's law. Their equa- 
tions (10), (13), (16) and (18) are the estimates of mo as a 
function of n and a number of known constants specific to 
the fits to observed aftershock sequences. We can use these 
relations of mo as a function of n to eliminate mo from equa- 
tion (12) to obtain direct estimates of n as a function of the 
measurable n a . For simplicity, we restrict the use of their 
findings to the case a = b. The estimate resulting from the 
fits performed by Helmstetter et al. [2004] yielded 

mo = TOmaB _(__)_ __ (14) 

with the values m m ax = 8.5, mj = 3, 6 = 0.1, c = 0.001, 
6=1 and Kf it = 0.008. The study initiated by Felzer et al. 
[2002] provided another estimate 
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where m max = 8.5, m d = 3, 6t = 0.08, At = 0.116 
days~ eT , b = a = 1, c = 0.014 and Mi = 6.04. Using 
the declustering performed by Reasenberg and Jones [1989], 
Sornette and Werner [2004] obtained 



mo = m ma x 



1-n 61n(10) 



(1 - 10 



-b(m max -m d )s 



(16) 



where m max = 8.5, m d = 3, 6 = 0.08, a = —1.67, c = 0.05 
and 6 = 1. Finally, using Bath's law, Sornette and Werner 
[2004] found 



mo = ra, 



( " ) 

V l-n ; 61n(10) 



-10 



b(Mx —m a ) 



(17) 



(15) 



where Mi — m a = 1.2 according to the the law, 6=1, 
m max = 8.5, and m d = 3. 

Substituting these four estimates of mo from equations 
(14), (15), (16), and (17) into equation (12) for n a provides 
four estimates of n a versus n all in terms of known constants. 
These four estimates of n as a function of n a can be used 
to find the correct fraction of aftershocks from the measur- 
able apparent fraction of aftershocks. Figure 4 shows these 
four estimates with the above constants. As noted above, n 
varies linearly with n a with a slope determined by mo. The 
four estimates of n as a function of n a allow for all possible 
values of mo. The solid line n = n a , corresponding to the 
slope 1 when mo = md, separates the left side of the graph, 
where mo < m d , from the right side, where mo > m d , which 
can be ruled out based on observed aftershocks from mag- 
nitudes 2. Thus, only the region to the left of the diagonal 
should be considered. 

Figure 4 can be used to find the real fraction of after- 
shocks from the measured apparent fraction by assuming 
one of the four estimates of mo as a function of n. For ex- 
ample, Helmstetter et al. [2004] find that 55 percent of all 
earthquakes are aftershocks above ma = 3. Using their val- 
ues to estimate mo as a function of n, we can determine that 
the real fraction of aftershocks is closer to 75 percent. Thus 
the size of this effect is significant. Furthermore, having de- 
termined a point on the line estimating n from n a for all 
values of mo fixes the slope of n(n a ) and therefore mo. Us- 
ing their values, we find mo = 1.2. Similar estimates can be 
made using the apparent fraction of aftershock values found 
by Felzer et al. [2002] and Reasenberg and Jones [1989]. 

Assuming that current maximum likelihood estimation 
methods of the ETAS model parameters, which assume 
mo = md, determine a branching ratio that corresponds 
to the present apparent branching ratio, we can similarly 
correct these values to find the true fraction of aftershocks 
using Figure 4. For example, Zhuang et al. [2004] find a 
"criticality parameter" of about 45 percent, which we take 
as a proxy for n a . Figure 4 shows that the true branching 
ratio then lies between 0.45 and 0.80, depending on which 
estimate (among the four models (14), (15), (16), and (17)) 
of mo as a function of n is chosen. These calculations sug- 
gest that previous estimates of the fraction of aftershocks 
obtained by various declustering methods significantly un- 
derestimated its value. 

3.2. Determination of apparent background events 
S a of uncorrelated seismicity 

In order to derive the number of shocks within one cas- 
cade that are not triggered by a mother above the thresh- 
old and thus appear as independent background events, we 
need to distinguish between the case where the initial (main) 
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shock of magnitude Mi is observable (i.e. Mi > m d ) and 
the case where it is undetected (i.e. Mi < m d ). 

If Mi > md, then the initial mainshock produces 
p(Mi)f bs observed direct aftershocks. On average, these 
will in turn collectively produce p(Mi) fob s n a observed sec- 
ond generation aftershocks. We specifically do not consider 
events above m d triggered from below m d , which we deal 
with below in the definition of the apparent background 
sources. By continuing this "above-water" cascade for all 
generations of aftershocks, we can calculate the number of 
triggered events that are in direct lineage above the thresh- 
old back to the mainshock as the infinite sum of terms of 
p(M\)f oos multiplied by the apparent branching ratio n a to 
the power of the generation. If, on the other hand, the initial 
mainshock is below raj, then no such direct "above- water" 
cascade will be seen. Any observed shock will be triggered 
by an event below the water. Thus, for the two cases, the 
"above-water" sequence is expressed as: 



bove 



P(Mi)U, 
1-na 



= N, 



obs i_ 







, Mi > m d 
,M 1 <m d 



(18) 



Furthermore, since in the ETAS model, a small earthquake 
may trigger large earthquakes, an event below m d may pro- 
duce an observed event above m d . An inversion method 
that reconstructs the entire branching structure of the model 
from an earthquake catalog will identify these shocks as 
background events. But since in reality these events were 
triggered by earthquakes below the detection threshold, we 
will refer to them as apparent background events. These 
events can of course trigger their own cascades. We thus 
define apparent background source S a as the number of ob- 
served events above m d that are apparently not triggered, 
i.e. have "mothers" below m d . Again, we distinguish be- 
tween the cases where the mainshock magnitude is Mi > m d 
and Mi < m d . For the first, S a is given by the total number 
of aftershocks below the threshold multiplied by the aver- 
age number r of direct aftershocks they trigger above the 
threshold. For the second case, we must also include the di- 
rect aftershocks of the initial mainshock that are observed: 



S a 



p(Mi) 



(1 - fobs)r ,Mi>m d 
H^ 1 - fobs)r + p(Mi)/ oi)s ,Mi<m d 



p(Ml) 



(19) 



Now, the number r of observable direct aftershocks above 
m d averaged over unobserved mothers between mo and m d 
is given by the following conditional branching ratio: 



m„ 



P(m\m < m d )p(m)f bsdm 



= (n - n a ) 



fob 



1- fob 



(20) 
(21) 



where we have used P(m\m < m d ) = P(m)/(1 — fobs) for 
m < m d and zero otherwise. Substituting (21) into the ex- 
pression for the apparent source (19) and re-arranging using 
(8), we obtain 



£ Y^-fobs{n - n a ) ,M 1 >m d 
bs , Mi < m d 
N bs(n - n a ) , Mi > rn d 



Sa 



4^/06^ -n a )+p(Mi)/ 



N bs(n — n a )+ p{M-i)f bs 



, Mi < m d 



(22) 



Equation (22) shows that, for each genuine background 
event, a perfect inversion method would count S a apparent 
background events. Figure 5 plots the number of apparent 
background events S a as a function of the branching ratio 



n for an example aftershock cascade set off by a magnitude 
m — 5 initial shock. We assumed the values m max = 8.5, 
m d = 3 and a = b = 1.0. The figure shows that for one 
cascade, i.e. one independent background event, hundreds 
of earthquakes appear as apparent background events when 
m < ma- 
in Figure 6, we investigate the relative importance of the 
apparent background events with respect to the observed 
number of aftershocks of one cascade. According to equa- 
tion (22) 



Sa/N bs = n-n a , 



(23) 



i.e. a significant fraction n — n a of events of the actually ob- 
served aftershocks are falsely identified as background events 
(since all events are really triggered from a single mainshock 
in our example). For m d — mo, the ratio is zero, since no 
events trigger below the detection threshold. However, as 
mo decreases and more and more events fall below m d , the 
fraction increases until n a goes to zero and the ratio ap- 
proaches n. This effect increases with decreasing a. Small 
values of a generally place more importance on the cumula- 
tive triggering of small earthquakes. 

3.3. Consistency check: N b s as the sum of "above- 
water" cascades triggered by the mainshock and by 
the apparent background events 

To complete the calculations and show consistency of the 
results, we now demonstrate that the observed cascades set 
off by the apparent background events, when added to the 
original "above- water" cascade, add up to the total observed 
number of aftershocks of the whole sequence. Each apparent 
source event will trigger its own cascade above the threshold 
m d with branching ratio n a ■ The total number of events due 
to the apparent background events and their cascades above 
the threshold is 



S a + S a Tl a + S a Tl a + .. 



S a 



1 — n a 

Substituting expression (22) and using (8) gives 



(24) 



N obl 



(n — n a J 



AT ("-"■) 1 p(M 1 -m )/ ob3 



, Mi > m d 
,Mi<m d 



(25) 



Combining the direct "above- water" cascade (18) with the 
apparent source cascades (25) gives the total amount of ap- 
parent events observed after the initial event 



N a 



■^V sou 



ce ~~f~ -t* above 
AT , ("-"a) 



p(M 1 -m )f obs 



(„- 



N o»- l-„, 



1) 1 p(.M 1 -m )f obs 1 q 
1-na 



Nobs, 



,M 1 >m d 
, Mi <m d 
(26) 



where N b s is given by (8). The last equality confirms the 
consistency of our decomposition into apparently-triggered 
earthquakes and apparent sources. 

Expressions (10) and (22) show that analyzing the tree 
structure of triggered seismicity only above the detection 
threshold leads to the introduction of an apparent source S a 
and an apparent branching ratio n a . It is important to real- 
ize that both are renormalized simultaneously by m d ^ mo- 
If the current inversion techniques for the ETAS parameters 
were perfect and correctly reconstructed the tree structure of 
all sequences, the inverted values would be equal to our an- 
alytical results (10) and (22). Accordingly, the value of the 
background source would be overestimated and the branch- 
ing ratio underestimated. In fact, one single true sequence 
will appear as many different sequences, each apparently set 
off by an apparent background event. 
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4. Conclusions 

We have shown that unbiased estimates of the fraction 
of aftershocks and the number of independent background 
events are simultaneously renormalized to apparent values 
when the smallest triggering earthquake mo is smaller than 
the detection threshold m^. In summary, mainshocks above 
the threshold will appear to have less aftershocks, resulting 
in a smaller apparent branching ratio. Meanwhile, unob- 
served events can trigger events above the threshold giv- 
ing rise to apparently independent background events that 
seem to increase the constant background rate to an ap- 
parent rate. Assuming that current techniques which are 
used to invert for the parameters of the ETAS model (for 
example, the maximum likelihood method) under the as- 
sumption md = mo are unbiased estimators of n a and S a , 
then the obtained values for the fraction of aftershocks and 
the background source rate correspond to renormalized val- 
ues because of the assumption that the detection threshold 
ra<j equals the smallest triggering earthquake mo. We pre- 
dict that n will be drastically underestimated and S strongly 
overestimated for mo much smaller than m d . 
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Figure 1. Schematic representations of the branching 
structure of the real ETAS model (left) and the apparent 
ETAS model (right). The initial mainshock is circled. 
Only events above the detection threshold rrid are ob- 
served. The apparent branching ratio does not take into 
account unobserved triggered events (dashed lines). An 
observed event triggered by a mother below mj appears 
as an untriggered background source event (circled). 




fraction of aftershocks n 

Figure 2. The apparent fraction of aftershocks (ap- 
parent branching ratio) n a varies linearly with the real 
fraction of aftershocks (real branching ratio) n with a 
slope fixed by the smallest triggering earthquake mo. As 
mo decreases, the apparent fraction of aftershocks signifi- 
cantly underestimates the real fraction. As examples, we 
chose mo = m d = 3 (solid), i.e. n a = n and no events 
are missed; m = (dashed); m = —5 (dotted); and 
mo = — 10 (dash-dotted). We further assumed parame- 
ters ma = 3, m m ax = 8, b = 1, and a = 1.0. A small 
value of a amplifies this effect (see Figure 3) . 
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Magnitude m of the smallest triggering earthquake 

Figure 3. The ratio of the apparent fraction of after- 
shocks (apparent branching ratio) n a over the real frac- 
tion of aftershocks (real branching ratio) n varies as a 
function of the smallest triggering earthquake mo. For 
mo = nid, n a — n and all events are detected above the 
threshold. For a small value of mo, the ratio becomes 
small, indicating that n a significantly underestimates n. 
Decreasing a amplifies this effect. We used parameters 
rrid = 3 (vertical reference line), 6=1. We 

varied a = 0.5 (dash-dotted), a = 0.8 (dashed), a — 1.0 
(solid). 
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Figure 4. The fraction of aftershocks (branching ra- 
tio) n can be estimated from the apparent fraction of 
aftershocks (apparent branching ratio) n a by using four 
estimates of the smallest triggering earthquake mo as 
a function of n as determined in Sornette and Werner 
[2004] (see text). The estimates of mo as a function of n 
were obtained from comparisons of the ETAS model pre- 
diction of the number of observed aftershocks and fits to 
observed aftershock sequences performed by Helmstetter 
et al. [2004] (solid), Felzer et al. [2002] (dash-dotted), 
Reasenberg and Jones [1989] (dotted) and from Bath's 
law (dashed). The additional diagonal solid line n a — n 
corresponds to mo = m d (no undetected events). Along 
any of the four lines, mo varies from minus infinity to 
rrimax- Given that we can rule out mo > ma, we can 
restrict the physical range to the left side of the refer- 
ence curve n a = n. Given an estimate of the apparent 
(measured) fraction of aftershocks n a , we can estimate 
the real fraction from one of the four lines. Because n is 
linearly proportional to n a , an estimate of n a not only 
determines n from one of the curves, but also estimates 
mo from the slope of the line connecting the point (n a , n) 
to the origin. For example, the estimate of Helmstetter et 
al. [2004] of 55 percent of observed aftershocks approx- 
imately gives a real fraction n of 75 percent according 
to their own fits to aftershocks and thereby determines 
mo = 1.2 for m ma x = 8 and md = 3. 
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Fraction of aftershocks n 

Figure 5. The number of apparent background events 
S a in an aftershock cascade due to a single background 
event of magnitude Mi = 5 as a function of the frac- 
tion of aftershocks (branching ratio) n for several values 
of the smallest triggering earthquake. For mo = m d , 
no events are missed. Therefore the number of appar- 
ent background events is zero. As mo decreases, events 
below the detection threshold trigger events above the 
threshold and hence the number of apparent background 
events increases. We vary mo = mj = 3 (solid, coinciding 
with x-axis), mo = (dash-dotted), mo = —5 (dashed), 
and mo = — 10 (upper solid curve). We used parame- 
ters m d = 3, nimax = 8, b = 1, and a — 1.0. For very 
small mo and n close to 1, almost all events above the 
detection threshold are triggered from below and thus S a 
becomes very large (see Figure 6). This effect is amplified 
for decreasing a (not shown). 
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Figure 6. The ratio of the number of apparent back- 
ground events S a over the total observed number Nobs of 
aftershocks of one cascade varies as n — n a . Here, we show 
the ratio as a function of the branching ratio n by assum- 
ing a particular value of mo. For mo = ma (solid, coin- 
ciding with x-axis), there are no apparent background 
sources. For mo less than mj, the ratio increases as more 
and more of the observed events are triggered by unob- 
served events. As examples, we show the ratio S a /N i, s 
for m = m d — 3 (solid, coinciding with x-axis), m = 
(upper solid line), mo = —5 (dashed) and mo = — 10 
(dash-dotted) as a function of the branching ratio n (av- 
erage number of aftershocks per earthquake also equal 
to the fraction of aftershocks in a catalog) for parame- 
ters md = 3, m max = 8, b = 1, and a = 1.0. For very 
small mo, n a approaches zero and the ratio S a /N b s ap- 
proaches its limiting value n, meaning that almost all 
observed earthquakes were triggered by events below the 
detection threshold m d . The effect of unobserved events 
triggering observed quakes resulting in an apparent back- 
ground source rate is further amplified by smaller values 
of a (not shown). 
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